Chaos and rectification of electromagnetic wave in a lateral semiconductor superlattice 
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We find the conditions for a rectification of electromagnetic wave in a lateral semiconductor 
superlattice with a high mobility of electrons. The rectification is assisted by a transition to a 
dissipative chaos at a very high mobility. We show that mechanism responsible for the rectification 
is a creation of warm electrons in the superlattice miniband caused by an interplay of the effects of 
nonlinearity and finite band width. 
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The research of nonlinear transport properties of the 
strongly coupled semiconductor superlattices (SLs) un- 
der the action of constant and alternating fields is the 
subject of traditional interest associated with the great 
potential of numerous applications in microwave and ter- 
ahertz technologies In recent years the primary goal 
has been investigations of novel strongly nonlinear phe- 
nomena at the miniband transport in SLs, such as dissi- 
pative chaos 1^ and excitation of Bloch oscillations under 
the action of a pure ac electric field ||, |. The last effect 
can be considered as a new mechanism of the rectifica- 
tion of THz radiation; the generated dc voltage is either 
integer or fractionally quantized [Q . 

The recent experimental progress in the observation of 
the influence of intensive THz electromagnetic radiation 
on voltage-current characteristics of SLs has stimu- 
lated advances in the theory ||, |[ ^. Observations of 
strongly nonlinear properties of the miniband transport 
in SLs subjected to an intensive electromagnetic field re- 
quire a relatively high level of the carrier density together 
with the high carrier mobility H, ^ ^. However, due to 
the experimental conditions associated with bulk SLs it is 
difficult to meet simultaneously these two requirements, 
because of a relatively low electron mobility and a for- 
mation of high field domains inside heavy doped SLs. 

These limitations can be overcome in the lateral semi- 
conductor superlattices (LSSLs) 1^, and especially in 
those LSSLs which are fabricated using the cleaved edge 
overgrowth (CEO) technique 0. In these LSSL devices a 
two-dimensional electron gas resides in an atomically pre- 
cise one-dimensional superlattice potential. CEO LSSLs 
are most suitable for an observation of Bloch oscillations: 
they have high electron mobility and can deliver high 
currents without formation of high field and low field do- 
mains 1^. In these respects the LSSLs are promising 
objects for investigations of various strongly nonlinear 
phenomena. 

In this Letter, we describe a novel mechanism of rectifi- 
cation of electromagnetic radiation in a superlattice that 
is different from the mechanism of AC-induced Bloch os- 
cillations ||] and that always results in a generation of an 
unquantizcd bias. We consider an interaction of a plane 



electromagnetic wave with 2D electrons in a single mini- 
band of LSSL, and incorporate the feedback effects of 
electron motion on the incident field. These feedback ef- 
fects are responsible for the nonlinearity of the problem. 
We show that for some field strengths and frequencies, 
the interplay of effects of the nonlinearity and the finite 
miniband width causes a weak asymmetry of a distribu- 
tion function of electrons in the momentum space with 
respect to the center of Brillouin zone. This deformation 
of nonequilibrium distribution function creates a direct 
current that is the reason of an appearance of bias in 
the wave transmitted through the LSSL. In mathemati- 
cal description, this effect corresponds to the symmetry- 
breaking bifurcation in nonlinear dynamical systems [|j . 
We also show that the interaction of the electromagnetic 
radiation with miniband electrons of the very high mo- 
bility can result in a dissipative chaos instead of a recti- 
fication. 

Consider a plane electromagnetic wave Ein = 
EQCOs(flt — kz), which is polarized along the superlat- 
tice axis, being incident normally to the surface of LSSL 
located at z = 0. The LSSL sits on a semi-insulating 
background with an average dielectric constant eo. We 
assume that both the superlattice thickness in the z- 
direction and the superlattice length along its axis are 
much less than the characteristic scale of the electromag- 
netic wave in the medium, 2Tr/k. To study the electron 
transport through a single miniband, a spatially homo- 
geneous LSSL with the period a and the miniband width 
A, we use the tight-binding energy-quasimomentum dis- 
persion relation e{kei) — (A/2) [1 — cos(fce/a)] (fee/ is the 
electron wave vector along the superlattice axis). The 
dynamics of miniband electrons is described by the fol- 
lowing nonlinear balance equations |^ 



V ~ —Uw — 7W, 
w — Uv — 'y{w — We 

U = Um(t) - KV, 



(la) 
(lb) 
(Ic) 



where J7i„(t) = {2ea/h{l + ))Ein{t,z = 0) and 
U{t) = {ea/h)Eout{t,z = 0) {Ein{t,z) and Eout{t, z) are 
the fields of the incident wave and the wave transmit- 
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ted by LSSL, respectively). The relaxation processes 
for miniband electrons are characterized by a scatter- 
ing constant 7. The electron variables v — moVa/h, 
w = {e — A/2)(A/2)~^ and Wf.q are the scaled elec- 
tron velocity, the scaled electron energy, and the equi- 
librium value of the scaled electron energy, respectively, 
and mo = (2fi.^)/(Aa^) is the effective mass at the bot- 
tom of the miniband. The lower (upper) edge of the 
miniband corresponds to w = —1 (w = +1), and the 
value of Weq {weq < 0) is a function of lattice tempera- 
ture. The scaled variables v{t) and w{t) are proportional 
to the variables V(t) and e{t), which are the electron ve- 
locity and the energy averaged over the time-dependent 
distribution function satisfying the Boltzmann equation. 

The first two Eqs. of set (|l|) are well-known super- 
lattice balance equations jlO|; their solution for a given 
field U (t) is exactly the same as the solution of the time- 
dependent Boltzmann equation for a tight-binding lat- 
tice. Eq. ( |la| ) describes the acceleration of electrons un- 
der the action of the electric field U{t) and their slowing 
down caused by an effective friction due to the scattering, 
while the second Eq. (lb) describes the balance of elec- 
tron's energy gain under the action of the electric field 
and the energy loss due to scattering. These equations 
are valid if the standard conditions of semiclassical ap- 
proach, eEoa < A, huj < Ag, eEoa < Ag {Ag is the 
width of minigap) 0|, are satisfied. 

The Eq. (|l^) is derived by modeling the LSSL as an 
equivalent current screen of infinitesimal thickness with 
the use of the Maxwell equations and appropriate bound- 
ary conditions |9| Such approach is valid if super- 
lattice thickness is much less than the wavelength 27r/fc 
1^; this condition is easily satisfied for typical LSSL ir- 
radiated by submillimeter or millimeter electromagnetic 
waves. This Eq. (0) describes the back influence of the 
electron current in LSSL on the dynamics of miniband 
electrons making the whole set of Eqs. (|l|) nonlinear. The 
degree of nonlinearity is controlled by the value of param- 
eter K — {An e'^Ns)/ 17100(1 + e^^), where c is the speed of 
light and Ns is the areal electron density. 

The dynamics of electrons in LSSL as well as the time- 
dependence of the transmitted and reflected waves is de- 
termined by an interplay of the relaxation (cx 7) and the 
nonlinearity (oc k Ng). Therefore, it is convenient to 

1 /2 

introduce a single parameter, F = (7/^) , controlling 

the relative contribution of these two basic processes. For 

1/2 

this we introduce the scaled time r = uJot, luq = [nj) , 
scaled fields u = U/uJo, Um = Um/oJo, and rewrite Eqs. 
(i as 

V = - (mi„(t) - F^-^u) w - Fu, 

W = {urn{T) -T^^v) V -T{w - Weq), (2) 

where the overdot means the differentiation with respect 
to the value r, Umir) ~ uqCoslut with uj = fl/tuo and 

1/2 

Mo = (2eai?o)/ftwo(l + ). The wave transmitted by 




FIG. 1: (color). The regions of the rectification (black) and 
the chaos (red) in the frequency - field strength plane for 
r = 0.1. Straight blue lines mark roots of Bessel function. 
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FIG. 2: (color). Same as in Fig.|l], but for F = 0.2. 



LSSL is u{t) = —T^^v + uq cos ujt. 

We have performed a systematic numerical search for 
stationary solutions {t ^ F~^) of Eqs. (||) which have 
a nonzero dc component of the outgoing wave, (u) 
(the averaging is made over a period of the ingoing wave 
T = 2Tr /uj). We also detect the chaotic attractors char- 
acterized by a positive value of the maximal Lyapunov 
exponent In contrast, the Lyapunov exponent is al- 
ways zero for the limit cycles. The locations of chaotic 
attractors and symmetry-broken periodic attractors with 
> 10"'' in the witQ-plane for different values of F are 
shown in Figs. |l[]|. The structure of the wwo-plot shown 
in Fig. |l] for F = 0.1 is quite typical for the behavior of 
our system at weak damping. It consists of the stripes 
for solutions supporting rectification and some area of 
chaotic solutions located at the lower values of uq. For 
a weak damping and a low frequency, "symmetry-broken 
stripes" are located near the lines of zero order Bessel 
function roots defined as Jo (uo /i^) = (Fig. |^, only 
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about 30 roots are shown). With an increase of damping 
(r = 0.2, FigJI), the region for an existence of chaotic 
attractors in the wMQ-plane shrinks, while the regions of 
nonchaotic attractors with (u) ^ occupy a few wide 
stripes (cf. Figs.|l| and |^). Finally, with a further in- 
crease of r, chaos completely disappears and attractors 
with (m) 7^ are located in a single compact region in 
the WMo-plane (FigJ|). For F > Tcr ~ 0.4 we didn't ob- 
serve the effect of rectification anymore. The generated 
dc voltage is unquantized (Fig.^. This is in contrast to 
the case of rectification making use of the mechanism of 
ac-induced Bloch oscillations, where the ratio (u) juj can 
be either an integer or a fractional number for a weak 
damping The efficiency of transformation of ac 

field to dc field, {u)/uq, is rather low; in conditions of 
Fig.^ it is less than four percents overall. 
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FIG. 3: The region of the rectification in the frequency - 
field strength plane for V — 0.3. 
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FIG. 4: The dependence of the spontaneously generated dc 
bias {u)/uj on frequency cu for F — 0.3. 

Physical mechanism responsible for the rectification of 
electromagnetic wave in LSSL will be our focus in the 
remainder of this paper. There are three main distinct 
physical factors which can contribute to the generation 
of unquantized dc voltage in our problem: finite band 
width, nonlinearity, and dissipation. In order to distin- 
guish an influence of these factors, we consider the limit 



F ^ in the Eqs. (||) and get 

^ -I- F^^ sin 6' = — uo coscjT, (3) 

where v = — sinO, w — — cos6 and u — ~6. Eq. (|^) de- 
scribes a time evolution of electrons within the first Bril- 
louin zone {\9\ < tt) for a negligible scattering but with 
an account of nonlinearity. Unbiased AC field can not 
create rotational states with {9) 7^ in the overdamped 
pendulum, therefore the AC-induced Bloch oscillations 
Q can not exist here. However, the ac field can excite 
symmetry-broken swinging states with 6*0 = (6*) 7^ [||. 
In a small vicinity of the symmetry-breaking bifurcation, 
a stationary solution of Eq. (||) can be presented in the 
form 9 — 00 + Q;sin(ct;T -|- /i), where a and /i are 
constants dependent on the parameters of LSSL. Sub- 
stituting this expression in Eq. (||) and equating zero 
harmonics, we have 

Jo(a)sin6lo = 0. (4) 

Eq. (U) shows that either 9o — {a. symmetric solution) 
or 9o ^ 0, but a should satisfy Jo (a) = (a boundary of 
symmetry-breaking, cf. ref. [p^). Analogously, equating 
first harmonics, we can find the dependence of amplitude 
a on uq and uj for a symmetric solution, as well as the 
dependence 9o on uq and uj at the boundary of symmetry- 
breaking ||l^. This analysis shows that the behavior of 
symmetric and symmetry-broken solutions is quite differ- 
ent: For a symmetric solution a increases with increase 
of Uq, while for an asymmetric regime 9o increases but a 
stays constant. 

Therefore, we should distinguish two regimes of the 
wave transmission through a LSSL. In the first regime, 
the interaction of ac field with an electron gas can led 
only to a difference between the amplitudes of incident, 
uo, and transmitted, Umax = ^^a, waves due to a non- 
linear dependence of a on uq {cf. The value of 6*0 
is zero indicating that the distribution function of elec- 
trons in fce/-space is symmetric in respect to the center 
of Brillouin zone. With an increase of uq this regime 
continues until a reaches the value corresponding to one 
of the Bessel roots. In this second regime of a strong 
interaction between electrons and radiation, a variation 
of amplitude of the incident wave practically doesn't re- 
sult in a change of amplitude of the transmitted wave: 
Umax oc a=const independently on uq, at least near the 
boundary of such regime. Instead, the energy of absorbed 
radiation goes to an effective heating of electron gas and 
a creation of "warm electrons" . The distribution func- 
tion of these "warm electrons" is slightly asymmetric in 
respect to the center of Brillouin zone (^o / 0). In the 
presence of scattering the asymmetry of distribution func- 
tion results in a directed current and in an appearance of 
dc component of the outgoing wave. 

Consider the physics of a rectification of the wave with 
a low frequency, w <^ 1, and a large amplitude, uq > F~^, 
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in more detail. In this case our analysis shows |13| that 
a — > uq/lu and therefore the condition for symmetry- 
breaking bifurcation in Eq. (||) becomes Jq (uq/lo) = 
(we also have checked it solving Eq. (^) numerically) . On 
other hand, this condition is in a reasonable agreement 
with the results of numerical simulations of rectification 
at small T employing balance equations (^) (see Figj^ for 
r = 0.1). Further, we find that the large amplitude of 
ac field is almost unchanged after passing trough LSSL, 
Umax = ss ui (uq / uj) = uq. The absolute value of gen- 
erated dc bias is very small: ~ 10"'* for F = 0.1 
(however, such values are still quite distinguishable nu- 
merically from artificial dc effects arising due to round 
off errors that result in < 10~*). The expression 
■^0 {uq/lj) — determines the condition for dynamic lo- 
calization of miniband electrons driven by a given ac 
field [|l^ . Both quantum mechanical [l5[| and semiclassi- 
cal [|16[ considerations of the dynamic localization show 
that a forward motion of an electron is transformed to 
its local oscillations. These localized electrons strongly 
interact with radiation. The localization of every elec- 
tron under the action of a strong ac field is a coherent 
control effect [Q. However, different oscillating elec- 
trons have different and random phases in the presence 
of scattering; interaction of many electrons with the ra- 
diation results in an effective electron's heating ||l^. It 
was shown recently ||l^ that the overheating may lead to 
the symmetry-breaking instability resulting in dc. The 
development of this instability and its next stabilization 
is a rather complicated process which is associated with 
a nonlinear feedback action of electrons on the incident 
field (the terms +T~'^vw and —T~^v^ in Eqs. (^). Im- 
portantly, the nonlinear effects can not be ignored any- 
more in the conditions of localization even for a low elec- 
tron density Ns- Thus, when a dissipation is weak and 
the self-consistent field inside LSSL is only slightly dif- 
ferent from the incident ac field, just an interplay of the 
effects of dynamic localization and nonlinearity provides 
a conversion of radiation into a very small dc bias. 

For available LSSLs fabricated with the use of CEO 
technique, a = 15 nm, A = 12 meV, Ag = 85 meV, 
eo = 13 (GaAs), iV^ 6 x 10" cm'^ and = 10 ps at 
temperature 0.3 K 0, we estimate F = 0.6 > Vcr ~ 0.4. 
Because F oc (A^^rAa^)"*/^, a further increase of the 
electron's density, the mobility, as well as an increase of 
the miniband width should result in a possibility to ob- 
serve the predicted effects of rectification and chaos. In 
physical units, the desirable field strength u'^^'^ « 2.75 
and the frequency oj « 0.35 (cf. Fig. |3|) correspond 
to e'^'^^ « 825 V/cm and f7/27r w 9 GHz, that can 
be easily reached within the range for existing Gunn 
diodes. At the same time, the conditions of applicabil- 
ity of the semiclassical superlattice equations are well 
satisfied for such values of parameters (in particular, 
ea£;n/A « uo(fiwo/A) ~ 10"^ < 1). 



In conclusion, we should notice that likely similar ef- 
fects of a resonant photovoltaic response of 2D electron 
gas in semiconductor hetero junctions to microwave and 
far-infrared radiations have been observed in the exper- 
iments |l|, |l|, 1^. These photovoltaic effects were at- 
tributed in Refs. |l^, ^ to a "resonant heating" 
of intraband electrons confined by the different methods, 
such as an excitation of the discrete 2D plasmons [ p^ , 
an application of the quantized magnetic field p9[ or a 
creation of the field-effect confined quantum dots ]2C|| . In 
contrast, we consider the localization of electrons in the 
superlattice miniband induced by the ac electric field it- 
self (dynamic localization). Possible realizations of the 
symmetry-breaking bifurcations in the particular exper- 
imental situations jl^, |l9|, |2^ must be further investi- 
gated. 
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